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We investigate the phase structure of lattice QCD for the general number of flavors in the param- 
eter space of gauge coupling constant and quark mass, employing the one-plaquette gauge action 
and the standard Wilson quark action. Performing a series of simulations for the number of flavors 
Nf = 6-360 with degenerate-mass quarks, we find that when Nf > 7 there is a line of a bulk first 
order phase transition between the confined phase and a deconfined phase at a finite current quark 
mass in the strong coupling region and the intermediate coupling region. The massless quark line 
exists only in the deconfined phase. 

Based on these numerical results in the strong coupling limit and in the intermediate coupling 
region, we propose the following phase structure, depending on the number of flavors whose masses 
are less than Ad which is the physical scale characterizing the phase transition in the weak coupling 
region: When Nf > 17, there is only a trivial IR fixed point and therefore the theory in the 
continuum limit is free. On the other hand, when 16 > Nf > 7, there is a non-trivial IR fixed 
point and therefore the theory is non-trivial with anomalous dimensions, however, without quark 
confinement. Theories which satisfy both quark confinement and spontaneous chiral symmetry 
breaking in the continuum limit exist only for Nf < 6. 

PACS numbers: 12.38.Gc 
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I. INTRODUCTION 

The fundamental properties of QCD are quark confine- 
ment, asymptotic freedom and spontaneous breakdown 
of chiral symmetry. Among them, the asymptotic free- 
dom is lost when the number of flavors exceeds 16^. Thus 
the question which naturally arises is what is the con- 
straint on the number of flavors for quark confinement 
and/or the spontaneous breakdown of chiral symmetry. 

As asymptotic freedom is the nature at short distances, 
one can apply the perturbation theory to investigate the 
critical number for it. However, because quark con- 
finement and spontaneous chiral symmetry breaking are 
due to non-perturbative effects, one has to apply a non- 
perturbative method throughout the investigation of the 
critical numbers for them. We employ lattice QCD for 
the investigation in this work, since lattice QCD is the 
only known theory of QCD which is constructed non- 
perturbatively. 

Lattice QCD is a theory with fundamental parameters, 
the gauge coupling constant g and quark masses m q , de- 
fined on a lattice with lattice spacing a. The inverse of 
the lattice spacing a" 1 plays a role of an UV cutoff. In 
order to investigate properties of the theory in the con- 
tinuum limit, one has to first clarify the phase structure 
of lattice QCD at zero temperature for general number 
of flavors, and then identify an UV fixed point and/or an 
IR fixed point of an RG (Renormalization Group) trans- 
formation. When the existence of such fixed points is 
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established, one is able to conclude what kind of theory 
exists in the continuum limit. 

In our previous work it was shown that, even in the 
strong coupling limit, when the number of flavors Np is 
greater than or equal to seven, quarks are deconfined 
and chiral symmetry is restored at zero temperature, if 
the quark mass is lighter than a critical value which is of 
order a -1 . 

In this work we extend the region of the gauge cou- 
pling constant to weaker ones, and investigate the phase 
diagram for general number of flavors Np at zero tem- 
perature. We mainly consider the case where Nf quark 
masses are degenerate. However, we also discuss the non- 
degenerate case. 

We employ the simplest form of the action in this work, 
that is, the Wilson quark action and the standard one- 
plaquette gauge action. At finite lattice spacings, the 
Wilson quark action is not chirally symmetric even at 
vanishing bare quark mass due to the Wilson term which 
is added to a naive discretized Dirac action in order to 
lift the doublers. Because of the fact that the action does 
not hold chiral symmetry, the phase diagram becomes 
complicated in general. 

As the phase diagram turns out to be complicated on 
the lattice, we first give a brief summary in terms of 
the f3 function in Sec. [H] Further, the phase diagram 
for general number of flavors when the theory would be 
chirally symmetric is conjectured in Sec. IIIII This con- 
jecture is based on our proposal for the phase diagram 
for the lattice action we employ, and is given because it 
may help the reader to understand the whole structure of 
the phase diagram. After showing the brief summary in 
this way, the action we employ is given with some basic 
notations in Sec. IIVI Here some important facts related 
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to chiral property of the quark are also discussed. Then 
the main part of the paper, our proposal for the phase 
structure for the general number of flavors in the case of 
the Wilson quark action, is summarized in Sec.[V] before 
giving the detailed numerical results which lead to our 
proposal. We also comment on previous results obtained 
with the staggered quark action. After giving numerical 
parameters for our simulations in Sec. IVII numerical re- 
sults in the strong coupling limit and at finite coupling 
constants are given in Sees. IV 111 and IV111I respectively. 
All the investigations up to this point are limited to the 
degenerate quark mass case. In Sec. IIXI we extend the 
study to the non-degenerate case to discuss implication 
for physics. We give a summary in Sec. |X] In an Ap- 
pendix, we study the case of N c = 2 in the strong cou- 
pling limit. Preliminary results of our study have been 
presented in Refs. UQ. 



II. PHASE STRUCTURE IN TERMS OF BETA 
FUNCTION 

The perturbative beta function of QCD with Np fla- 
vors of massless quarks is universal up to two-loop order: 
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The coefficient bo changes its sign at Np — 16^. Hence, 
for Np > 17, asymptotic freedom is lost and the point 
g = is an IR fixed point. The theory governed by 
this IR fixed point is a free theory and the quark is not 
confined. The second coefficient b\ changes its sign at 
Nf ps 8.05. This implies, if one would take the two- 
loop form of the beta-function, that an IR fixed point 
appears at finite coupling constant for 9 < Np < 16. Of 
course we cannot trust the two-loop form at finite cou- 
pling constants. However, since the IR fixed point from 
the two-loop beta function locates in a weak coupling re- 
gion for Np ~ 16, it is plausible that the beta function 
has a non-trivial IR fixed point for N' < Np < 16 with 
some N' < 16. When such an IR fixed point exists, the 
coupling constant cannot become arbitrarily large in the 
IR region. This implies that quarks are not confined. 
The long distance behavior of the theory is governed by 
the non-trivial IR fixed point: It is a theory with an 
anomalous dimension. 

A pioneering study on the Np dependence of the QCD 
vacuum was made by Banks and Zaks in 1982 4]. Based 
on the result of the quark confinement in a pure gauge 
theory they assumed that the quark is confined and 




FIG. 1: Renormalization group beta function, (a) Conjecture 
by Banks and Zaks |4J assuming confinement in the strong 
coupling limit for all Nf- (b) Our conjecture deduced from 
the results of lattice simulations. 



that the beta function is negative in the strong coupling 
limit for all Np- Using the perturbative results men- 
tioned above, they conjectured Fig. ^a) as the simplest 
Np dependence of the beta function, and studied the 
phase structure of QCD based on this Np dependence of 
the beta function. Because of an additional non-trivial 
UV fixed point for Np > N', their conjecture for the 
phase structure is complicated. 

In the argument of Banks and Zaks, the assumption 
of confinement and negative beta function in the strong 
coupling limit plays an essential role. However, there 
exist no proofs of confinement in QCD for general Np 
even in the strong coupling limit. A non-perturbative 
investigation on the lattice is required. Therefore, in our 
previous work we performed numerical simulations 
of QCD in the strong coupling limit for various Np, us- 
ing the Wilson fermion formalism for lattice quarks. We 
found that, when Np > 7, quarks are deconfined and 
chiral symmetry is restored at zero temperature even in 
the strong coupling limit, when the quark mass is lighter 
than a critical value. 

In some previous literatures, quark confinement in the 
strong coupling limit has been assumed explicitly or im- 
plicitly. However, the arguments for the confinement in 
the strong coupling limit are based on either the large 
N c limit |3, a meanfield approximation (1/d expansion) 
0, , or a heavy quark mass expansion 0, EJ ■ Because 
we are interested in the theory with dynamical quarks 
with N c kept 3, and effects of dynamical quarks become 
significant only when they are light, the results of these 
approximations cannot be applied. Actually, our numer- 
ical result at g = oo explicitly shows violation of quark 
confinement at small quark masses when Np > 7. 

Here we extend the study to weaker couplings. Based 
on the numerical results obtained on the lattice (see the 
following sections for details) combined with the results 
of the perturbation theory, we conjecture Fig. ^b) for 
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the Np dependence of the beta function: When Np < 6, 
the beta function is negative for all values of g. Quarks 
are confined and the chiral symmetry is spontaneously 
broken at zero temperature. (Corresponding critical Np 
is 2 for the case N c = 2 See the discussion in the 
Appendix.) On the other hand, when Nf is equal or 
larger than 17, we conjecture that the beta function is 
positive for all g, in contrast to the conjecture by Banks 
and Zaks shown in Fig.^a). The theory is trivial in this 
case. When Np is between 7 and 16, the beta function 
changes sign from negative to positive with increasing g. 
Therefore, the theory has a non-trivial IR fixed point. 

Here we note that there are several related works which 
do not use lattice formulation These 
works gave interesting results for the quark confinement 
condition. However, the critical number for the confine- 
ment is at most suggestive, because we need fully non- 
perturbative method to investigate the problem. That 
is, the theory should be constructed from the beginning 
non-perturbatively. In this regards, lattice QCD is the 
only known theory constructed non-perturbatively. 



III. CONJECTURE FOR THE PHASE 
DIAGRAM IN CHIRALLY SYMMETRIC CASE 

In this section, we present our conjecture for the phase 
diagram of QCD, for the case when the theory is chiral 
symmetric. The conjecture is based on the phase diagram 
we propose for the case of the Wilson fermion action, 
where the chiral symmetry is violated. The reason why 
we make this kind of conjecture is that it may help the 
reader to understand the phase diagram of the Wilson 
fermion action which is more complicated. The phase 
diagrams we conjecture in the g — m,Q plane are presented 
in Fig. H 



A. N F < 6 

The phase diagram is simple. At zero temperature all 
the region is in the confined phase. See Fig. |2I a )- 
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B. 



N F > 17 



At zero temperature, there is a first-order phase tran- 
sition which separates a deconfined phase from the con- 
fined phase. The transition occurs at xn'^ %t%cal which is 
of the same order of magnitude as the inverse of a typ- 
ical correlation length £ of gluon dynamics, such as the 
plaquette-plaquette correlation length. 

In the confined phase, the nature of the system is es- 
sentially same as that for Np < 6: When mo ^> 
dynamical effects by quark loops can be neglected. In 
this case, the system is equivalent to the quenched QCD. 
That is, for mo 3> the phase is the confined phase. 
When dynamical effects of quark loops change the phase 



confinement 



N f <l6 



a) S 




g = 
IR 



de confinement 








first order phase transition 




confinement 





N F > 17 



b) 8 




g = 
UV 




deconfinement 
first order phase transition 
confinement 



7<:N F <\6 



c) s-- 



g = 



FIG. 2: Phase diagram for a chirally symmetric case, (a) 
Nf < 6, (b) N F > 17, (c) 7 < N F < 16. 



from that in the quenched case, there will be a phase 
transition in general, and this kind of the transition oc- 
curs at mo ~ £ _1 . 

When the gauge coupling decreases towards zero, the 
dimension- less correlation length £/<z increases and di- 
verges at zero coupling constant in the confined phase. 
Therefore, along the phase transition line, the dimension- 
less critical quark mass should behave as m < £ %t%cal a = 
Ada. Here a vanishes towards g = 0, and A^ ~ £ _1 is the 
physical scale which characterizes the phase transition in 
the weak coupling region. 

The "?" mark in the weak coupling region means 
that we have performed numerical simulation for only the 
strong and intermediate coupling regions because of tech- 
nical reason. In this sense, our proposal for the strong 
and intermediate coupling regions is based on our nu- 
merical results, while our conjecture in the weak region 
is based on the assumption that the transition occurs at 

^critical t — 1 

In the upper region above the phase transition line 
quarks are not confined. There is only an IR fixed point 
at g = 0. That is, the theory in the continuum limit is 
free. 
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We would like to make a comment that there is 
an alternative possibility that the transition occurs at 
a -1 , where a -1 is the inverse of the lat- 



m 



critical 



tice spacing. In the strong and intermediate coupling 
regions is of order a -1 . Therefore, there is no es- 
sential difference between the two possibilities whether 
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critical 




In the usual argument of the decoupling theorem for 
QCD, the effect of the particle whose mass is much heav- 
ier than the QCD scale parameter Aqc d can be absorbed 
by renormalization of physical quantities. Based on a 
similar consideration, one may argue that the particles 
whose masses are much heavier than Aqc p> are irrelevant 
for quark confinement. This is the reason why we assume 
that the transition occurs at TO g"** ca/ ^ However, 
whether a theory exists in a sense of constructive field 
theory is a different problem. Therefore, we think that it 
is not possible to disregard the alternative possibility of 



^critical 



~ l only from a theoretical argument. 



Nonetheless we think that the assumption that the 



transition occurs at to 



critical 




is more plausible than 
the assumption m <£ iUcal ^ a^ 1 . In the following we as- 
sume that the phase transition occurs at TO g"* lca ' ~ 
This is our conjecture. 

However, we also refer to the alternative possibility of 
^critical ^ a -i w \ ien j s appropriate to refer to. In 
this alternative case, the transition point at g = is off 
the massless point and therefore all finite massive quarks 
belong to the deconfined phase. 



It will be possible to see in future which of m ^ %tlcal 



a 



or to 



critical 



£ 1 is realized by numerical simula- 



tions. 



C. 7 < N F < 16 

The phase diagram at zero temperature shown in 
Fig. |5fc) is similar to that in the case of Np > 17. The 
phase transition line locates around mo ~ due to 
the same reason as in the case of Np > 17. 

However, the flow of RG transformation is different. 
The g = point in this case is an UV fixed point, con- 
trary to the case Np > 17. On the other hand, quarks 
arc deconfined in the upper region above the phase tran- 
sition line as in the case of Np > 17. Therefore there 
should be an IR fixed point at finite coupling constant. 
Otherwise gauge coupling would become arbitrary large 
at large distance, and quarks would be confined in the 
strong coupling limit contrary to the result of the nu- 
merical simulation. 



IV. FUNDAMENTALS OF LATTICE QCD AND 
THE ACTION 



Lattice QCD is defined on a hypercubic lattice in 4- 
dimensional euclidean space with lattice spacing a. A 
site is denoted by a vector n — (ni, 112, n^, n^j, where 



rii's are integers. A link with end points at the sites n 
and n + jl is specified by a pair (n, //), where ft denotes a 
unit vector in the /i direction. 



Action 



The gauge variable C/„ iAl which is an element of SU(3) 
gauge group is defined on the link (n, fi). The action for 
gluons which we adopt is given by 

Sgauge = -J ^] ^ r (^n,^^n+A,^^n+/i+i>,/j ' (^) 



9 



where g is the gauge coupling constant. This action is 
called the standard one-plaquette gauge action. We usu- 
ally use, instead of the bare gauge coupling constant g, 
(3 defined by 
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(5) 



The quark variable is a Grassmann number defined on 
each site. It is well-known that a naive discretization of 
the Dirac action 



Sfe 



QnQn i 



(0) 



with too being the bare fermion mass, leads to 16 poles 
instead of one pole. This problem is called the species 
doubling. To avoid this problem K. Wilson proposed to 
add a dimension 5 operator called the Wilson term 



'^2(q n q n +p, - 2q n q n + qn+p,q n ), 



(7) 



to the naive discretized action. Rescaling the fermion 
field by q — \/2Kip and making the action gauge invari- 
ant, we obtain the Wilson fermion action 

Swilson = a 3 2^[VViVn ~ K {^n(l ~ T/J^n.juVVi+A 
n 

+^ n+lx (l+ 7fi )UlM}. (8) 



where 



K = 



1 



2(TO a + 4) 



(9) 



which is called the hopping parameter. 

The full action S is given by the sum of the gauge part 



~>gauge and the fermion part S 

Wilson i 



/=i 



(10) 
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where / is for flavors. 

The expectation value of an operator 0(U,ip,-ip) is 
given by 

n,fj, nj 

(11) 

where Z is the partition function 

Z = J II dU ^ II d^ /) d^/)exp(5), (12) 
with dUn^ being the Haar measure of SU(3). 

B. Fundamental parameters 

In the case of degenerate Np flavors, lattice QCD 
contains two parameters: the gauge coupling constant 
(3 = 6/g 2 and the bare quark mass or the hopping pa- 
rameter K = l/(moa+4). In the non-degenerate case, we 
have, in general, Np independent bare masses (hopping 
parameters). 

In this work we consider mainly the case where Np 
quark masses are degenerate, because it is simpler. How- 
ever, the conjecture can be extended to non-degenerate 
cases. We will discuss this point in Sec. IIXI 



C. Continuum limit 

As mentioned in the Introduction, in order to see what 
kind of theory in the continuum limit exists, one has 
to investigate the phase structure and identify UV fixed 
points and IR fixed points. 

When the g — point is an UV fixed point, the theory 
governed by this fixed point is an asymptotically free 
theory. On the other hand, when the g = point is an 
IR fixed point, the theory governed by this fixed point 
is a free theory. When there is a IR fixed point at finite 
coupling constant g, the theory is a non-trivial theory 
with the long distance behavior determined by the IR 
fixed point. 



D. Finite temperatures 

At finite temperatures the linear extension in the time 
direction N t is much smaller than those in the spatial 
directions (N x ,N y ,N z ). The temperature T is given by 
1/N t a. For gluons the periodic boundary condition is im- 
posed, while for quarks an anti-periodic boundary condi- 
tion is imposed in the time direction. 

Although we are interested in the phase structure of 
lattice QCD at zero temperature, we investigate phase 
transitions on lattices at finite temperatures. Increasing 
the value of Nt , we carefully examine the Nt dependence 



of the phase transition. If the transition is N t indepen- 
dent for sufficiently large Nt , then the transition is a bulk 
transition (phase transition at zero temperature). 



E. Quark mass and chiral symmetry 

In the Wilson quark formalism, the flavor symmetry 
as well as C, P and T symmetries are exactly satisfied 
on a lattice with a finite lattice spacing. However, chiral 
symmetry is explicitly broken by the Wilson term even 
for the vanishing bare quark mass too = at finite lat- 
tice spacings. The lack of chiral symmetry causes much 
conceptual and technical difficulties in numerical simula- 
tions and physics interpretation of data. (See for more 
details Ref. [l4[ and references cited there.) 

The chiral property of the Wilson fermion action was 
first systematically investigated through Ward-Takahashi 
identities by Bochicchio et al. 15]. We also indepen- 
dently proposed to define the current quark mass 
by 

2m,(0|P|7r) = -m 7r (0| A 4 \ir), (13) 

where P is the pseudoscalar density and A4 the fourth 
component of the local axial vector current. We use this 
definition of the current quark mass as the quark mass 
in this work. In general we need multiplicative normal- 
ization factors for the pseudoscalar density and the axial 
current. In this study, we absorb these renormalization 
factors into the definition of the quark mass, because 
this definition is sufficient for later use. We note that the 
quark mass thus defined has an additive renormalization 
constant to the bare quark mass, because the Wilson 
quark action does not hold chiral symmetry. 

With this definition of quark mass, when the quark is 
confined and the chiral symmetry is spontaneously bro- 
ken, the pion mass vanishes in the chiral limit where the 
quark mass vanishes at zero temperature. 

However, in the deconfined phase at finite tempera- 
tures the pion mass does not vanish in the chiral limit. It 
is almost equal to twice the lowest Matsubara frequency 
n/Nt. This implies that the pion state is approximately 
a free two-quark state. The pion mass is nearly equal to 
the scalar meson mass, and the rho meson mass to the 
axial vector meson mass; they are all nearly equal to the 
twice the lowest Matsubara frequency ir/N t . Thus the 
chiral symmetry is also manifest within corrections due 
to finite lattice spacing. 

In a previous study [T3 we found that, given the gauge 
coupling constant and bare quark masses, the value of the 
quark mass does not depend on whether the system is in 
the high or the low temperature phase when the gauge 
coupling is not so strong: (3 > 5.5 for Np — 2. However, 
this property is not guaranteed for smaller [3. In general, 
the quark mass in the deconfined phase differs from that 
in the confined phase in the strong coupling region. 
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FIG. 3: (a) m q at j3 = oo. (b) m-n at ft = oo. Results with 
an anti-periodic boundary condition (apbc) in the t-direction 
and those with the periodic boundary condition (pbc) are 
compared on Nt — 8 and 4 lattices. 



0.25 



N F <6 



K 




0.125 



confined 



m.q = oo 



P 



0.25 



K 



N F <6 



confined 




b) 



finite temperature 
QCD phase transition] 



P 



niq = oo 



0.125 



FIG. 4: The phase structure for Nf < 6; (a) at zero tem- 
perature, and (b) at finite temperatures. The chiral limit 
(massless quark limit) is shown by thick curves labeled by 
"m 9 = 0", and the finite temperature QCD transition at a 
fixed finite Nt is shown by a shaded curve. 
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F. Quark mass at g = 

As mentioned earlier, the Wilson term lifts the clou- 
biers and retains only one pole around mo = in the free 
case. On the other hand, there are other poles at differ- 
ent values of the bare masses. They are remnants of the 
doublers. 

In Fig. Eta) the quark mass defined by is plotted 
versus the bare quark mass for the case of free Wilson 
quark. At mo > (1/K > 8) the quark mass behaves 
as expected: it monotonously increases with m . On the 
other hand, at mo < (1/K < 8) the behavior of the 
quark mass is complicated: it does not monotonously 
decrease with decreasing mo, but increases after some 
decrease and becomes zero at a finite negative value of 
mo . Usually, the region of negative bare quark mass is ir- 
relevant for numerical calculations of physical quantities. 
However, this region is also important for understanding 
the phase diagram. 

We also plot in Fig. EJb) the mass of the pion which 
is composed of two free quarks with periodic boundary 
condition for the time direction and with anti-periodic 
boundary condition. These results will be referred to 



V. PHASE DIAGRAM 

Here we propose, based on our numerical results which 
will be shown later and the perturbation theory, the 
phase diagram in the (j3, K) plane for the Wilson quark 
action coupled with the one-plaquette gauge action. 



A. iV F < 6 

For Np < 6, the chiral line K c ((3) where the cur- 
rent quark m q vanishes, exists in the confined phase(See 
Fig. Eta)). The value of K c ((3) at (3 = oo is 1/8, which 
corresponds to the vanishing bare quark mass mo = 0. 
As (3 is decreased, the value of K c increases up to 1 /4 at 
j3 = 0. If the action would be chiral symmetric, the chiral 
line should be a constant, 1/8, as in Sec. IIIII The line 
K = corresponds to the case of infinitely heavy quarks. 
Quarks are confined for any value of the current quark 
mass for all values of (3 at zero temperature (N t = oo). 

On a lattice with a fixed finite Nt, we have the finite 
temperature deconfming transition at finite /3, because 
the temperature T = 1/N t a becomes larger as (3 in- 
creases in asymptotically free theories. At K = (m q = 
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FIG. 5: Phase diagram for Nf = 240. Dark shaded (green) 
lines represent our conjecture for the bulk transition line in 
the limit Nt — oo. Points connected by dashed lines are the 
measured points where quarks are massless. Light shaded 
lines are for the finite temperature transition at Nt — 4 and 



oo), the first order finite temperature phase transition 
of pure SU(3) gauge theory locates at (3 C = 5.69254(24) 
and 5.89405(51) for N t = 4 and 6 18] and at /3 C ~ 6.0625 
for Nt — 8 19]. This finite temperature transition turns 
into a crossover transition at intermediate values of K, 
and becomes stronger again towards the chiral limit K c . 
As K is increased, the finite temperature transition line 
crosses the K c line at finite (3 uM- We note that, for 
understanding the whole phase structure which includes 
the region above the K c line (negative values of the bare 
quark mass), the existence of the Aoki phase is impor- 
tant [2C|. A schematic diagram of the phase structure 
for this case is shown in Fig. 0£b). F° r simplicity, we 
omit the phase structure above the K c line. It is known 
that the system is not singular on the K c line in the 
high-temperature phase (to the right of the finite tem- 
perature transition line) [14J . The location of the finite 
temperature transition line moves toward larger j3 as Nt 
is increased. In the limit Nt — oo, the finite tempera- 
ture transition line will shift to [3 — oo so that only the 
confined phase is realized at T = 0. 



B. When Nf is very large 

We present the result for the case of Np — 240 in 
Fig. [SJ The reason why we investigate the case where the 
number of flavor is so large as 240 is the following: We 
have first investigated the case of Np = 18 as a generic 
case for Nn > 17. However it has turned out that the 
phase diagram looks complicated when Np — 18. So, to 
understand the phase structure for Np > 17, we have 
increased the number of flavors like 18, 60, 120, 180, 
240, and 300, and systematically viewed the results of 
the quark mass and the pion mass for all these numbers 
of flavors. Then we have found that when the number of 
flavors is very large as 240, the phase diagram is simple as 
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b) 

FIG. 6: The same as Fig. but for (a) N F = 18 and (b) 
N F = 12. 



the chirally symmetric case discussed in Sec. IIIII There- 
fore we first show the result for the case of Nf = 240. 

At finite Nt where numerical simulations have been 
performed, the finite temperature transition occurs as 
shown in Fig. [S] As Nt increases, the transition line 
moves towards larger value of (3. The envelop of those 
finite temperature transition lines is the zero temperature 
phase transition line, shown by the dark shaded curve in 
the figure. Note that, at finite Nt, the system is not 
singular on the part of the dark shaded line in the right 
hand side of the finite temperature transition line (i.e., 
in the high temperature phase). 

The salient fact is that, at zero temperature, the mass- 
less line exists only in the deconfined phase and passes 
through from f3 = down to /3 = oo, i.e., it is quite sim- 
ilar to the chirally symmetric case of Nf > 17 shown in 
Fig. OJb) . Thus the IR fixed point at g = governs the 
long distance behavior of the theory and therefore the 
theory is a free theory. 



C. Nf > 17, but not so large 

A typical phase diagram for Np > 17 looks like that in 
Fig. Eta), where the case of Nf = 18 is displayed. The 
massless line in the deconfined phase which starts from 
(3 = oo hits the first-order phase transition around p = 



8 



4, and goes underground crossing the first-order phase 
transition line. 

The place where the crossing occurs moves towards 
(3 = as the number of flavors increases, and finally 
reaches the (3 = axis when Np > 240. The massless 
line exists only in the deconfined phase and it starts from 
(3 = oo (g = 0). Therefore the IR fixed point at g = 
governs the long distance behavior of the theory. That 
is, the theory is free. 

D. 7 < N F < 16 

The phase diagram for Np = 12 is shown in Fig.[|J{b), 
which looks similar to that of Np = 18. The salient fact 
is that the massless line exists only in the deconfined 
phase. That is, there is no massless line in the confined 
phase. Therefore, in the continuum limit, the quark is 
not confined. 

The difference from the case of Np > 17 is that the 
g = point is an UV fixed point in this case. Thus there 
should be an IR fixed point at finite coupling constant. 
If there would be no IR fixed point, we would encounter 
a contradiction that on one side quarks are not confined, 
but on the other side the gauge coupling constant be- 
comes arbitrary large as the distance between quarks be- 
comes large. 

When we increase the number of flavors from Np = 
6 to 17, the form of the (3 function changes from that 
in the upper frame (Np < 6) in Fig. ^b) to that in 
the lower frame (Np < 17) through that in the middle 
frame (7 < Np < 16). We safely assume that the form 
changes smoothly for varying Np. That is, for Np = 7 
the IR fixed point should appear at very large coupling, 
and gradually the position of the IR fixed point moves 
towards the g = point. For Np = 16, the IR fixed point 
is expected to be close to g = 0. 

We have been unable to identify numerically the posi- 
tion of the IR fixed point. Technically it is not easy to do 
so. It might be even one of the metastable states beyond 
the first-order phase transition line, i.e., under the first 
sheet of the phase diagram. 



TABLE I: Simulation parameters for Nf = 6 on lattices 8 2 x 
10 x N t (N t = 4, 6, 8), 12 3 x N t (N t = 12) and 18 2 x 24 x N t 



(N t 


= 18). 






P 


K 


At 


4 


0.0 


0.2, 0.21, 0.22, 0.235, 0.24, 0.25 


n ni-n nn95 




0.1 


0.2495 


n ni 




0.2 


0.249, 0.24936 


n ni 

U.U -L 




0.3 


0.2485, 0.249 


01 

U.U -L 




0.4 


0.248 


01 

U.U -L 




O.o 


n oo n ooc n o/i n O/ic n 0/17C 
U.zo, U.zoo, U.24, 0.24o, 0.24ro 


n 01 

U.U -L 




i n 
1.0 


no n oi n oo n ooc: n oo n ook 
0.2, U.21, 0.22, 0.220, 0.2o, U.2oO, 








0.2of, 0.2oo, U.24, U.24o 


01 -0 005 




o < \ 
Z.U 


n o/i 
U.24 


0.01 




i n 

4.U 


n oo 
U.22 


0.01 




4.0 


U.lO, U.lD, U.lDO, U.lDD, U.lD/, U.lDO, 








U.loy, U.l(, U.lo, U.ly, U.214o 


0.01 


r 
u 


0.0 


0.2, 0.21, 0.22 


01 

U.U -L 




O.O 


0.24(^0 


n 01 

U.U -L 




1.0 


n a c 
0.240 


0.01 




2.0 


n /i 
0.24 


0.01 




4.0 


n 00 
0.22 


0.01 




4.0 


n it; nice nic? n ico ni 1 ? n 1 71 
O.lD, U.lDO, U.lD/, U.lDO, U.l/, U.l/1, 








n 1 71 e n 1 ti t n 1 71 n n 1 70 n 1 70c 
O.l/lO, 0.1/1/, 0.1/ly, 0.1/2, 0.1/20, 








a 1 70 r\17/1 A 1 7r nio 

0.173, 0.174, 0.175, 0.18 


0.01 


Q 
O 


0.5 


0.2475 


01 

U.Ul 




i n 
l.U 


U.z4o 


0.01 




z.U 


U.Z4 


0.01 




4.U 


o 00 


0.01 




1 EC 

4.0 


nit; n i n i n 17 n 1 70 n 1 70Q 
U.lD, U.lO/, U.lDO, U.l<, U.liz, U.lizo, 








U.l<24, U.l/25, U.l/3, 0.1/0, O.lo 


0.01 




5.5 


0.1615 


0.01 


12 


4.5 


0.175 


0.01 


18 


0.3 


0.2485 


0.01 




0.4 


0.248 


0.01 




0.5 


0.2475 


0.01 




1.0 


0.245 


0.01 




4.5 


0.172, 0.1725, 0.173, 0.175, 0.2143 


0.01 



E. Previous studies using staggered quarks 

We note that there are several works investigating the 
critical number of flavors for quark confinement using 
the staggered quarks. The staggered fermion (Kogut- 
Susskind fermion) (2^ is a formulation of lattice fermions 
different from the Wilson fermion we adopt. Unlike the 
Wilson fermion action ©, the staggered fermion action 
explicitly violates the flavor symmetry due to flavor- 
mixing interactions at finite lattice spacings. However, 
because a part of flavor-chiral symmetry is preserved on 
the lattice, a numerical analysis of chiral properties is 
easier than the Wilson fermion. On the other hand, the 
staggered fermion action can describe quarks only when 



Np is a multiple of 4. A trick to study the cases Np 7^ 4n 
is to modify by hand the power of the fermionic determi- 
nant in the numerical path-integration. This necessarily 
makes the action non-local which sometimes poses con- 
ceptually and technically difficult problems. 

Evidences of strong first order transition are reported 
with staggered quarks for N F = 4-18 HI iH S 
I27I l2£j | . Chiral symmetry is restored at the transition. 
Furthermore, the transition was shown to be a bulk phase 
transition for several cases of Np = 8 0123, 12 H and 
16 [2j3. 

The most systematic study was done by the Columbia 
group in Ref. [27|. They studied the case Np = 8 at a 
fixed bare quark mass mod = 0.015 on 16 3 x Nt lattices 
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TABLE II: The same as Table Q] but for N F = 7. 


Nt 


P 


K 


At 


4 


0.0 


0.2, 0.21, 0.22, 0.23, 0.235, 0.24, 0.245, 0.25, 0.2857, 0.3333 


0.01-0.005 




1.0 


0.245 


0.01 




2.0 


0.125, 0.1429, 0.1667, 0.2, 0.2083, 0.2174, 0.2273, 0.24, 0.25 


0.01 




3.0 


0.235 


0.01 




4.0 


0.125, 0.1429, 0.1667, 0.178, 0.185, 0.195, 0.2, 0.2226, 0.25 


0.01 




4.5 


0.14, 0.15, 0.16, 0.161, 0.162, 0.163, 0.164, 0.165, 0.17, 0.18, 0.19, 0.2143, 0.25 


0.01 




5.0 


0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.1982, 0.21 


0.01 




5.5 


0.1, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.161 


0.01 




6.0 


0.08, 0.09, 0.1, 0.11, 0.1111, 0.12, 0.125, 0.13, 0.14, 0.1429, 0.15, 0.1519, 0.152, 








0.155, 0.16, 0.1667, 0.1724, 0.1786, 0.2, 0.25, 0.3333 


0.01 


6 


0.0 


0.2, 0.21, 0.22, 0.23, 0.24, 0.245, 0.25 


0.01 




2.0 


0.193, 0.214, 0.24 


0.01 




4.0 


0.163, 0.178, 0.195, 0.217 


0.01 




4.5 


0.15, 0.16, 0.164, 0.165, 0.168, 0.1681, 0.1682, 0.169, 0.17, 0.18 


0.01 




5.0 


0.12, 0.13, 0.14, 0.15, 0.16, 0.17 


0.01 




5.5 


0.12, 0.13, 0.135, 0.138, 0.14, 0.15 


0.01 


8 


0.0 


0.25, 0.2857, 0.3333 


0.01 




2.0 


0.1667, 0.2, 0.2083, 0.2174, 0.2273, 0.24, 0.25 


0.01 




4.5 


0.15, 0.16, 0.164, 0.165, 0.168, 0.1683, 0.1684, 0.169, 0.18, 0.25 


0.01 




5.0 


0.12, 0.13, 0.14, 0.15, 0.16, 0.17 


0.01 




5.5 


0.1, 0.11, 0.12, 0.13, 0.135, 0.138, 0.14, 0.145, 0.15, 0.16 


0.01 




6.0 


0.1, 0.11, 0.12, 0.13, 0.14, 0.145, 0.1476, 0.1519, 0.1667, 0.1724, 0.1786, 0.2, 0.25 


0.01 


18 


0.0 


0.245, 0.25 


0.01-0.005 




4.5 


0.165, 0.1675, 0.1684, 0.17, 0.18, 0.19, 0.2143 


0.01 




5.5 


0.135, 0.15, 0.162 


0.01 




6.0 


0.1, 0.11, 0.12, 0.13, 0.14, 0.145, 0.1476, 0.1519 


0.01 





TABLE III: 


The same 


as Table □ but for N F = 8, 10, 14, 17, 120, 180, and 360. 






Nt 


P 


K 


Ar 


8 


4 


0.0 


0.22, 0.23, 0.24, 0.25 


0.01 




18 


0.0 


0.25 


0.01 


10 


4 


0.0 


0.25 


0.01 


14 


4 


0.0 


0.25 


0.01 


17 


4 


6.0 


0.145, 0.1473, 0.15 


0.01 




8 


6.0 


0.1473 


0.01 




18 


6.0 


0.1473, 0.15 


0.01 


120 


4 


0.0 


0.125, 0.1316, 0.1389, 0.1429, 0.1471 


0.005-0.0025 


180 


4 


0.0 


0.122, 0.125, 0.1274, 0.1303, 0.1333, 0.1429 


0.005-0.0025 


360 


4 


0.0 


0.125 


0.00125 




8 


0.0 


0.125 


0.00125 



where Nt = 4-32. They found that the transition point 
shifts towards larger (3 as Nt is increased from 4 to 8, 
but stays around C = 4.73 for N t — 8 and 16. They 
concluded that this is a bulk transition, and proposed an 
interpretation that this transition is an outgrowth of the 
crossover transition between the strong and weak cou- 



pling regions observed in pure SU(3) gauge theory. 

Although their interpretation is in clear contrast with 
our proposal we note that their numerical results them- 
selves are consistent with our phase diagram shown in 
Fig- El As an illustration, let us suppose that we perform 
simulations at, say, K = 0.14 in Fig. E{b). For small 
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TABLE IV: The same as Table Q but for N F = 12. 





P 




At 


4 


0.0 


0.18, 0.2, 0.205, 0.21, 0.215, 0.22, 0.23, 0.24, 0.25 


0.01 




2.0 


0.15, 0.17, 0.185, 0.19, 0.195, 0.2, 0.21, 0.25 


0.01 




4.0 


0.14, 0.15, 0.155, 0.16, 0.165, 0.17, 0.18, 0.2, 0.25 


0.01 




4.5 


0.13, 0.14, 0.145, 0.15, 0.16, 0.1667, 0.17, 0.1786, 0.1852, 0.2, 0.25 


0.01 




5.0 


0.11, 0.12, 0.13, 0.135, 0.14, 0.15, 0.16, 0.1667, 0.1724, 0.1786, 0.2, 0.25 


0.01 




6.0 


0.1, 0.1111, 0.125, 0.1429, 0.1667, 0.1724, 0.1786, 0.2, 0.25, 0.3333 


0.01 




10.0 


0.13, 0.14, 0.15, 0.16 


0.01 


6 


0.0 


0.18, 0.2, 0.21, 0.215, 0.22, 0.23, 0.24, 0.25 


0.01 




2.0 


0.15, 0.17, 0.185, 0.19, 0.195, 0.2, 0.21, 0.2156, 0.2253, 0.236 


0.01 




4.0 


0.14, 0.145, 0.15, 0.155, 0.16, 0.165, 0.17, 0.18, 0.1949, 0.2028, 0.2114 


0.01 




4.5 


0.13, 0.14, 0.145, 0.15, 0.155, 0.16, 0.17 


0.01 




5.0 


0.11, 0.12, 0.13, 0.135, 0.14, 0.145, 0.15, 0.16 


0.01 


8 


0.0 


0.2, 0.205, 0.2083, 0.21, 0.215, 0.22, 0.2222, 0.225, 0.23, 0.24, 0.25 


0.01-0.0025 




2.0 


0.15, 0.17, 0.185, 0.19, 0.195, 0.2, 0.21, 0.225, 0.25 


0.01 




4.5 


0.13, 0.14, 0.15, 0.155, 0.16, 0.1667, 0.1786, 0.2, 0.25 


0.01 




6.0 


0.1429, 0.1667, 0.2, 0.25 


0.01 


18 


4.5 


0.14, 0.15, 0.155, 0.16 


0.01 



TABLE V: The same as Table [fl but for N F = 16. 

N t §_ k At 

4 0.0 0.18, 0.19, 0.2, 0.21, 0.22, 0.25 0.01 
4.5 0.125, 0.135, 14, 0.1429, 0.145, 0.15, 

0.1563, 0.1667,0.2, 0.25 0.01 

6.0 0.125, 0.1429, 0.1493, 0.1667, 0.2, 0.25 0.01 

10.0 0.1111, 0.125, 0.1429, 0.1667, 0.2, 0.25 0.005-0.0025 

100. 0.1111, 0.125, 0.1429, 0.1667, 0.2, 0.25 0.005-0.0025 

8 0.0 0.18, 0.19, 0.2, 0.21, 0.22, 0.25, 0.27 0.01 
4.5 0.125, 0.1429, 0.145, 0.15, 0.155, 

0.1563, 0.1667, 0.2, 0.25 0.01 

6.0 0.125, 0.1429, 0.1493, 0.1667, 0.2, 0.25 0.01 

10.0 0.1111, 0.125, 0.1429, 0.1667, 0.2, 0.25 0.005 

100. 0.1111, 0.125, 0.1429, 0.1667, 0.2, 0.25 0.005 



18 0.0 


0.25 


0.01 


4.5 


0.15 


0.01 



Nt we have the finite temperature transition at finite 0. 
This transition point moves towards larger as we in- 
crease Nt, but, eventually, stops on the dark shaded line 
when Nt becomes larger than a critical value. In order 
to clarify the whole phase structure in the case of stag- 
gered quarks and to discriminate different possibilities, 
it is indispensable to perform a more systematic study 
exploring a wider region of the parameter space of quark 
mass and gauge coupling constant. 



VI. PARAMETERS FOR NUMERICAL 
SIMULATIONS 

We perform simulations on lattices 8 2 x 10 x Nt (N t = 4, 
6 or 8), 16 2 x 24 x N t (N t = 16) and 18 2 x 24 x N t 
(Nt = 18). We vary Np from 2 to 360, selecting some 
typical values of Np. For each Np, we study the phase 
structure in the coupling parameter space (0, K). Sim- 
ulations for Np < 6 are discussed in [lj, )2!m • We sum- 
marize simulation parameters for Np > 6 in Tables ITllXI 
where we list values of Nt, 0, K and At. (For readers 
who are interested in more details we will provide them 
on request.) We adopt an anti-periodic boundary condi- 
tion for quarks in the t direction and periodic boundary 
conditions otherwise. In the cases where a thermalized 
state is achieved, typical statistics for hadronic measure- 
ments are about 10-100 configurations sampled every 1— 
5 trajectories, where the length of one trajectory is one 
molecular-dynamics time. When the hadron spectrum 
is calculated, the lattice is duplicated in the direction of 
lattice size 10 for Nt < 8, which we call the z direction. 
Statistical errors are estimated by the jack-knife method. 

We use the hybrid R algorithm [2lJ for the genera- 
tion of gauge configurations. The R algorithm has dis- 
cretization errors of 0(NpAr 2 ) for the step size At of a 
molecular dynamic evolution. As Np increases we have 
to decrease At, such as At =0.0025 for N F = 240, to 
reduce the errors. We have checked that the errors in the 
physical observables we study are sufficiently small with 
our choices of Ar for typical cases. 

It should be noted that, in QCD with dynamical 
quarks, there are no order parameters for quark confine- 
ment. We discuss confinement by measuring the screen- 
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TABLE VI: The same as Table Q but for N F = 18. 


N t 




13 K 


At 


4 




0.0 0.17, 0.175, 0.18, 0.19, 0.195, 0.2, 0.205, 0.21, 0.215, 0.22, 0.235, 0.245, 0.25 


0.02-0.005 






2.0 0.16, 0.17, 0.18, 0.19, 0.2, 0.21, 0.25 


0.01 






4.0 0.13, 0.14, 0.15, 0.152, 0.154, 0.16, 0.1613, 0.1667, 0.17, 0.18, 0.2, 0.25 


0.01 






4.5 0.115, 0.125, 0.135, 0.14, 145, 0.15, 0.1613, 0.1667, 0.1724, 0.2, 0.2143, 0.25 


0.01 






6.0 0.1, 0.1111, 0.125, 0.1429, 0.1613, 0.1667, 0.1724, 0.2, 0.25 


0.01 






10.0 0.13, 0.14, 0.15 


0.01 


8 




0.0 0.18, 0.19, 0.2, 0.21, 0.22, 0.23, 0.25, 0.27 


0.01 






2.0 0.16, U.l(, 0.18, 0.19, U.2, 0.21 


U.01 






4.0 0.13, 0.14, 0.15, 0.152, 0.154, 0.156, 0.158, 0.16, 0.17, 0.18 


0.01 






4.5 0.135, 0.14, 0.145, 0.15, 0.155, 0.1667, 0.2, 0.25 


0.01 






6.0 0.125, 0.1389, 0.1429, 0.1667, 0.2, 0.25 


0.01 


18 




0.0 0.25 


0.01 






4.5 0.15, 0.2143 


0.01 


TABLE VII: The same as Table [fl but for N F = 60. 




P 


K 


At 


4 


0.0 


0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1429, 0.1538, 0.1613, 0.1667, 0.1724, 0.2, 0.25, 0.3333, 0.5, 1.0 0.01-0.00125 




6.0 


0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1333, 0.1429, 0.1493, 0.1538, 0.16, 0.1667, 0.2, 0.25 


0.005 


8 


6.0 


0.125, 0.1429, 0.1667, 0.2, 0.25 


0.005 


TABLE VIII: The same as Table Q but for N F = 240. For N t = 16 the spatial lattice is 16 2 x 24. 


Nt 




P K 


At 


4 




0.0 0.08333, 0.09091, 0.1, 0.1111, 0.122, 0.1234, 0.125, 0.1333, 0.137, 0.1429, 0.1493, 0.1538, 








0.1667, 0.2, 0.25 


0.0025 






2.0 0.09091, 0.1, 0.1111, 0.125, 0.1359, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025 






4.5 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1343, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025 






6.0 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1337, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025-0.00125 






100. 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1293, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025-0.00125 




1000. 0.125 


0.0025 


8 




0.0 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.129, 0.1333, 0.137, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025 






4.5 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1316, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025 






6.0 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1309, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025 






100. 0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1264, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.0025 


16 




0.0 0.125, 0.129, 0.1333, 0.1429, 0.1538, 0.2 


0.002-0.00125 


TABLE IX: The same as Tabled but for N F = 300. For Nt = 16 the spatial lattice is 16 2 x 24. 


Nt 


P 




At 


4 


0.0 


0.08333, 0.09091, 0.1, 0.1111, 0.1143, 0.1176, 0.1212, 0.125, 0.1429, 0.1493, 0.1538, 0.1667, 0.2, 0.25 


0.00125 




4.5 


0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1429, 0.1538, 0.1667, 0.2, 0.25 


0.00125 




6.0 


0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1429, 0.1493, 0.1538, 0.1667, 0.2, 0.25 


0.00125 




100. 


0.08333, 0.09091, 0.1, 0.1111, 0.125, 0.1429, 0.1493, 0.1538, 0.1667, 0.2, 0.25 


0.00125 


8 


0.0 


0.08333, 0.09091, 0.1, 0.1111, 0.1143, 0.1176, 0.1212, 0.125, 0.1429, 0.1667, 0.2, 0.25 


0.0025-0.00125 
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FIG. 7: Plaquette (a) and Polyakov loop (b) at j3 — as 
a function of 1/K for various Nf on Nt — 4 lattices. Lines 
are to guide the eyes. In this figure, approximately vertical 
lines leaning to the right mean that we have two state signals 
there. 



ing pion mass, the screening quark mass, the values of 
the plaquette and the Polyakov loop. See Sec. IVIII for 
details. In the following, we call the pion screening mass 
simply the pion mass, and similarly for the quark mass. 

The numerical calculations were performed on var- 
ious computers including the dedicated parallel com- 
puter QCDPAX and Fujitsu VPP500 at the University 
of Tsukuba, and HITAC S820/80 at KEK. 



VII. RESULTS OF NUMERICAL SIMULATIONS 

AT g = oo 

We first study the strong coupling limit f3 = (g = oo). 
In a previous work [lj we have shown the following. When 
Np < 6, we have only one confined phase from the heavy 
quark limit K = up to the chiral limit K c = 0.25. 
On the other hand, for Np > 7, we found a strong first 
order transition at K = < K c . We confirmed that 
this transition Kd is a bulk phase transition (transition 
at zero-temperature) by increasing the value of Nt up to 
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FIG. 8: The plaquette at /3 — as a function of 1/K for (a) 
Nf = 240 and (b) N F = 300, at Nt = 4, 8 and 16. Lines are 
to guide the eyes. 



N t = 18 for N F = 7. When quarks are heavy (K < K d ), 
both the plaquette and the Polyakov loop are small, and 
satisfies the PCAC relation oc m q . Therefore 
quarks are confined and the chiral symmetry is sponta- 
neously broken in this phase. We found that m q in the 
confined phase is non-zero at the transition point Kd, 
i.e., the chiral limit does not belong to this phase. On 
the other hand, when quarks are light (K > Kd), the 
plaquette and the Polyakov loop are large. In this phase, 
m w remains large in the chiral limit and is almost equal 
to twice the lowest Matsubara frequency ir/N t . This im- 
plies that the pion state is an almost free two-quark state 
and, therefore, quarks are not confined in this phase. The 
pion mass is nearly equal to the scalar meson mass, and 
the rho meson mass to the axial vector meson mass. The 
chiral symmetry is also manifest within corrections due 
to finite lattice spacing. 

In this paper, we extend the study to larger Nf- In 
Fig. [7| we show the results of the plaquette and Polyakov 
loop at Nt = 4 for Np = 7 - 300. Clear first order transi- 
tions can be seen at 1/K larger than 1/0.25 = 4. We then 
study the N t dependence of the results, as shown in Fig.[S] 
for the cases Np = 240 and 300. We find that, although 
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FIG. 9: Results of ml and 2m q at f3 = for 7V F = 18-300. 
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FIG. 10: The transition point 1/Kd at /? = versus Nf for 
N t = 4 and Nt > 8. For clarity, data at Nt = 8 for N F = 300 
is slightly shifted to a larger Nf in the figure. 



the transition point shows a slight shift to smaller 1/K 
when we increase Nt from 4 to 8, the transition stays at 
the same point for N t > 8. For N F = 240, 1/K d ~ 8.1(1) 
at JVt = 4 and l/JC d = 7.8(2) at N t = 8 and 16. A similar 
result was reported for Nf — 7 in 1]. We conclude that 
the transition is a bulk transition for Np > 7. 

Figure El shows the results of and 2m q , in the lat- 
tice units, at (3 = for various numbers of flavors. Wc 
clearly see that at the exactly same hopping parameter 
K = Kd where the plaquette makes a gap, the pion mass 
and the quark mass also make gaps. When the quark 
is heavy [K < Kd), the PCAC relation (mj oc m q ) is 
well satisfied. On the other hand, when the quark mass 
is smaller than the critical value, the 1/K dependence of 
the pion mass squared and the quark mass m q looks 
strange at first sight. However, when one compares this 
dependence of m q with that of the free quark case shown 
in Fig.|3fa), one easily notices that the 1/K dependence 
is essentially the same as that of the free quark. The 
quark mass vanishes at 1/K ~ 7-8 for Np — 60, 240, 




FIG. 11: Results of (a) plaquette and (b) Polyakov loop for 
various Nf obtained at f3 = 4.5 on Nt — 4 lattices. 



and 300. This corresponds to the fact that the free quark 
mass vanishes at 1/K = 8. The 1/K dependence of 
is also essentially the same as that of the free quark case 
shown in Fig. OJb) . We stress that the chiral limit where 
the quark mass vanishes does not exist in the confined 
phase. 

We show a map of the phase transition point K = Kd 
at (3 = in Fig. [TD| 



VIII. RESULTS OF NUMERICAL 
SIMULATIONS AT FINITE g 

We now extend the study to finite (3 (g < oo). Results 
of the plaquette and Polyakov loops at (3 = 4.5 and 6.0 
are plotted in Figs. 1111 and 1121 respectively, for various 
Np. We note that at f3 = 4.5, both the plaquette and 
Polyakov loop shows singular behavior at some quark 
mass for Np < 18, in contrast to the cases of Np > 240. 

In order to understand the structure of the phase dia- 
gram for each number of flavors, from now on we discuss 
each Np separately: We first intensively investigate the 
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FIG. 12: The same as Fig. El but for /3 = 6.0. 



cases Np = 240 and 300, and then decrease Np, because 
we found that the phase structure is quite simple for Np 
> 240. 



A. N F = 240 and 300 

Figure ITTA shows the results of and 2m q for Np = 
240 and 300, obtained at various values of (3 on the Nt — 
4 lattice. A very striking fact is that the shape of 
and 2m q as a function of 1 /K only slightly changes in the 
deconfincd phase, when the value of f3 decreases from oo 
down 0. Only the position of the local minimum of 
at 1/K ~ 8, which corresponds to the vanishing point of 
m, 
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FIG. 13: 7r4 and 2m 9 for (a) N F = 240 and (b) 300, obtained 
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slightly shifts toward smaller 1/Jf. The results for FIG. 14: The same as Fig-EJ but for A F = 240 and A t = ; 



A_p = 240 and 300 are essentially the same, except for 
very small shifts of the transition point and the minimum 
point of mj. We obtain similar results also for N t — 8 
(see Fig. 1141 . That is, the massless line in the deconfined 
phase runs through from (3 = 00 to [3 = 0. Thus we 
obtain the phase diagram shown in Fig. [S] 

From the perturbation theory, the m q = point at 



(3 = 00 is a trivial IR fixed point for Np > 17. The 
phase diagram shown in Fig. [S] suggests that there are 
no other fixed points on the m q — line at finite (3. 
In order to confirm this, we investigate the direction of 
the renormalization group (RG) flow along the m q — 
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line for Np = 240, using a Monte Carlo renormalization 
group (MCRG) method. 

We make a block transformation for a change of scale 
factor 2, and estimate the quantity A/3 = /3(2a) — (3(a): 
We generate configurations on an 8 4 lattice on the m q — 
points at (3 = and 6.0 and make twice blockings. We 
also generate configurations on a 4 4 lattice and make once 
a blocking. Then we calculate A/3 by matching the value 
of the plaquette at each step. From a matching of our 
data, we obtain A/3 ~ 6.5 at (3 = and 10.5 at (3 = 
6.0. The value obtained from the two-loop perturbation 
theory is A/3 ~ 8.8 at /3 = 6.0. The signs are the same 
and the magnitudes are comparable. 

It is known for the pure SU(3) gauge theory that one 
has to make a more careful analysis using several types 
of Wilson loop with many blocking steps to extract a 
precise value of A/3. We reserve elaboration of this point 
and a fine tuning of 1/K at each (3 for future works. For 
Np = 240, because the velocity of the RG flow is large, 
we will be able to obtain the sign and an approximate 
value of A/3 by a simple matching. 

This result implies that the direction of the RG flow 
on the m q — line at (3 = and 6.0 is the same as that 



at /3 = oo. This further suggests that there are no fixed 
points at finite (3. All of the above imply that the theory 
is trivial for N F = 240. 



B. 240 > N F > 17 

Now we decrease Np from 240. As discussed in 
Sec. I VIII the deconfined phase transition point K& de- 
creases with decreasing Np in the strong coupling limit 
(3 = 0. However, except for this shift of the bulk transi- 
tion point, the 1/K dependence of and m q are quite 
similar in the deconfined phase when we vary Np from 
300 down to 17, as shown in Fig-IHI The results at (3 = 6.0 
and 4.5 shown in Fig. 1151 indicate that the 1/K depen- 
dence of and m q are almost identical to each other, 
except for a small shift toward smaller 1/K as Np is 
decreased. These facts imply that the structures of the 
deconfined phase are essentially identical from Np = 17 
to 300. 

As we show below, a closer examination of the data 
shows that the massless quark line in the deconfined 
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FIG. 17: N F = 18: Results of ml and 2m, versus 1/K. (a) 
/8 = 4.0, (b) (3 = 4.5, and (c) (3 = 6.0. 



phase hits the phase transition line at finite (3 when Np 
is not so large as Nf < 60, while it runs through from 
(3 = oo to [3 — when iVp is very large such as 240 and 
300. 

For Nf — 18 we make simulations at (3 = 0.0, 2.0, 4.0, 
4.5, and 6.0 on N t = 4 and 8 lattices, as listed in Table lVll 
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FIG. 18: N F = 7: Results of ml and 2m, versus 1/if. (a) 
(3 = 0.0, (b) (3 = 2.0. 



Results for and to 9 are shown in Figs. H*6l and IT7I We 

note that, at (3 = 4.5, the and m q show characteristic 
behavior of massless quarks around 1/K — 6.5 in the 
deconfined phase, both on Nt = 4 and 8 lattices. Such 
massless point exists also at (3 = 6.0, but is absent at 
/3 = 0.0 and 2.0. 

At (3 = 4.0 we carry out detailed simulations around 
the transition point. On the N t = A lattice, the first order 
deconfining phase transition locates around K = 0.150 
(1/K — 6.67), and the massless quark point exists at 
K ~ 0.158 (1/K ~ 6.33) in the deconfined phase. When 
we increase Nt to 8, the phase transition point shifts to 
K = 0.154-0.156 (1/K = 6.49-6.41) where we observe 
two-state signals lasting more than 450 trajectories. We 
confirm that the points K = 0.152 and 0.158 belong 
to the confined and deconfined phase, respectively. In 
the deconfined phase the massless quark point exists at 
K ~ 0.155 (1/K ~ 6.45). See Fig. Efta). From this 
we conclude that the massless quark line in the decon- 
fined phase hits the first order bulk phase transition line 
around (3 = 4.0. The phase structure for Np — 18 is 
summarized in Fig.^a). 
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FIG. 20: N F = 12: Results of ml and 2m q versus 1/K. (a) 
P = 2.0, (b) P = 4.5. 



C. 16 > N F > 7 

As stressed several times, the quark confinement is lost 
for Np >7at/3 = 0(g = oo). We now intensively 
simulate the cases Np = 7, 12 and 16 at several finite 
values of (3. The simulation parameters for Np = 7, 12 
and 16 are listed in Tables ITT1 llVl and M 

In Figs. 1181 and 1191 and 2m q at selected values of 
(3 are plotted for Np — 7. Results for Np = 12 and 16 
are similar, as shown in Figs. 1201 and !21l 

For Np = 7, 12 and 16, we find that the 1/K depen- 
dence of and 2m q in the deconfined phase are similar 
to those shown in Fig. ^| for Np > 17 at all the val- 
ues of p. That is, they are essentially the same as those 
of a free quark state shown in Fig. Careful looking 
at the values of and 2m q at (3 — 4.5 for both cases 
reveals that the massless line hits the phase transition 
line around (3 = 4.5, at 1/K ~ 6.1 for Np = 7 and at 
1/K 6.5 for Np = 12. 

Together with the data of plaquette and Polyakov loop 
(see Fig. 1221 as a typical example), we obtain the phase 
diagram shown in Fig. EJb) for Np = 12. The phase 



diagrams for Np — 7 and 16 are similar. We find that 
the gross feature of these phase diagrams is quite similar 
to the case Np = 18 shown in Fig. Eta). 

We note again that when Np < 16, the point g = 
is an UV fixed point. This is in clear difference with the 
case of N F > 17. 



D. N F < 6 

As already reported in ref . Jl| , numerical results show 
that, at j3 = 0, the confined phase extends to the chiral 
limit. Simulations at finite values of s how that the 
phase structure is similar for Np — 2-6 [14|. (See Tabled 
for run parameters for Np = 6.) Thus we have the phase 
diagram shown in Fig. 0] 



IX. IMPLICATION FOR PHYSICS 

We now extend the discussion to the cases of non- 
degenerate quarks. In the case of degenerate masses, the 
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global structure of the phase diagram is the same at large 
P for Np > 7. As far as the quark mass is below a critical 
value, the system is in the deconfmed phase. Therefore, 
we expect that our proposal for the relation of the phase 
structure and the number of flavors remain unchanged 
for non-degenerate cases, when we redefine Nf as the 
number of flavors which satisfies the condition too < A^: 
If there exist more than seven quarks whose masses are 
lighter than A<$, quarks should not be confined. In na- 
ture quarks are confined, therefore the number of flavors 
whose masses are less than A^ should be equal or less 
than six. The scale A^ is numerically calculable and will 
be obtained in the future. 

This conclusion is based on our assumption that the de- 
confining phase transition occurs at m < ^ %t%caX ~ Now 
we make a comment on the alternative possibility that 
the deconfining phase transition occurs at mg"** ~ 
a -1 . If this would be correct, the number of flavors 
should be equal or less than six irrespective of their 
masses. Since six species of quarks have been discovered, 
this implies that there are no more species of quarks. 
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X. CONCLUSIONS 

We have investigated numerically the phase structure 
of QCD for the general number of flavors. Performing a 
series of simulations for degenerate quark mass cases em- 
ploying the one-plaquette gauge action and the standard 
Wilson quark action, we have found that when Np > 7 
there is a line of a first order phase transition between the 
confined phase and a deconfmed phase at a finite current 
quark mass in the strong coupling region and the inter- 
mediate coupling region. The massless quark line exists 
only in the deconfmed phase. 

Based on these numerical results in the strong coupling 
limit and in the intermediate coupling region, together 
with an assumption that the phase transition occurs at 
m cnticai ^ £— l j n £ ne wea ]j coupling region, where £ is a 
typical correlation length of gluons, we propose the fol- 
lowing phase structure. The phase structure crucially 
depends on the number of flavors, Np, whose masses are 
less than A^ ~ which is the physical scale which char- 
acterizes the phase transition: When Np > 17, there 
is only a trivial IR fixed point and therefore the the- 
ory in the continuum limit is free. On the other hand, 
when 16 > Np > 7, there is a non-trivial IR fixed point 
and therefore the theory is non-trivial with anomalous di- 
mensions, however, without quark confinement. Theories 
which satisfy both quark confinement and spontaneous 
chiral symmetry breaking in the continuum limit exist 
only for Np < 6. If there exist more than seven quarks 
whose masses are lighter than Ad, quarks should not be 
confined. In nature quarks are confined, therefore the 
number of flavors whose masses are less than A^ should 
be equal or less than six. The scale A^ is numerically 
calculable and will be obtained in the future. 

We have discussed that, only from a theoretical argu- 
ment, we cannot exclude an alternative possibility that 
the deconfining phase transition occurs at m ^ nttcal ^ 
a -1 . This corresponds to the case A^ = oo. If this would 
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FIG. 23: Results of physical quantities at g — oo and K — K c 
for SU(2) QCD. (a) the plaquette and Polyakov loop, and (b) 
m% and 2m q . 



TABLE X: Simulation parameters for SU(2) QCD at /3 = 0.0 
performed on 8 2 x 10 x N t lattices with N t = 4 and 8. The 
number of configurations, sampled every 5 trajectories, are 
10-22, except for the Nf = 2 runs at K — 0.25 which are not 
thermalized. 
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that the SU(2) beta function has a characteristics similar 
to the SU(3) one: In the case of SU(2), the asymptotic 
freedom is lost when Np > 11 and the 2-loop result for 
the critical Np is 5, instead of 17 and 8, respectively, 
for SU(3) Qj. The smaller numbers for critical Np are 
natural because the confining force is weaker in SU(2). 

Adopting the standard plaquette gauge action and the 
Wilson quark action, we have performed a series of simu- 
lations in the strong coupling limit following the strategy 
of the SU(3) case. Our simulation parameters are com- 
piled in Table 13 



be correct, the total number of flavors should be equal or 
less than six irrespective of their masses: Since six species 
of quarks have been discovered, this implies that there are 
no more species of quarks. This conclusions is based on 
an additional assumption that QCD only is responsible 
to quark confinement. If some other interactions would 
also affect the quark confinement at some scale A new , the 
number of quarks, whose masses are smaller than A new , 
should be six, in this case. 

Which of m ^ iUcal a~ x or m ^ rltlcal ^ jg realized 
can be investigated by numerical simulations in future. 
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Appendix: SU(2) QCD 

As an extension of the color SU(3) QCD case, we study 
color SU(2) QCD in the strong coupling limit. We note 



First, performing a simulation for Np — 8 and Nt — 4, 
we confirm that the deconfinmg transition occurs at 
K = 0.2-0.21 and that the chiral limit K = K c = 0.25 
is in the deconfincd phase. Decreasing Np gradually at 
K = K c , we find that the chiral limit remains in the 
deconfined phase down to Nf = 3. When we further 
decrease Np to 2, the number of the inversion Ni nv in 
CG iterations shows a rapid increase with molecular- 
dynamics time and finally exceeds 10,000 in clear con- 
trast with small numbers of O(10 2 ) for Np > 3. We 
conclude that the chiral limit is in the confined phase for 
Np = 2. Fig. I23f a) shows our results of physical quanti- 
ties versus Np at K = K c . To study the Nt dependence 
of the transition, we simulate on an Nt = 8 lattice for the 
critical case Np = 3. The stability of the result, shown 
in Fig. l23f bh confirms that the deconfining transition we 
observe is a bulk phase transition. 

At finite coupling constants we may expect results sim- 
ilar to those of SU(3). Thus, we conjecture in paral- 
lel to the SU(3) case that, when Np > 11, the theory 
is free. On the other hand, when 10 > Np > 3, the 
theory is non-trivial with anomalous dimensions without 
quark confinement. Theories which satisfies both quark 
confinement and spontaneous chiral symmetry breaking 
exists only for Np < 2. 
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